Faraday patterns in coupled one-dimensional dipolar condensates 
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We study Faraday patterns in quasi-one-dimensional dipolar Bose-Einstein condensates with para- 
metrically driven dipolar interactions. We show that in the presence of a roton minimum in the 
excitation spectrum, the emergent Faraday waves differ substantially in two- and one-dimensional 
geometries, providing a clear example of the key role of confinement dimensionality in dipolar gases. 
Moreover, Faraday patterns constitute an excellent tool to study non-local effects in polar gases, 
as we illustrate for two parallel quasi-one-dimensional dipolar condensates. Non-local interactions 
between the condensates give rise to an excitation spectrum characterized by symmetric and anti- 
symmetric modes, even in the absence of hopping. We show that this feature, absent in non-dipolar 
gases, results in a critical driving frequency at which a marked transition occurs between correlated 
and anti-correlated Faraday patterns in the two condensates. Interestingly, at this critical frequency, 
the emergent Faraday pattern stems from a spontaneous symmetry breaking mechanism. 



PACS numbers: 03.75.Kk, 89.75.Kd, 05.30.Jp 

I. INTRODUCTION 

Inter-particle interactions play an essential role in the 
physics of ultra-cold gases. Although in many experi- 
ments these interactions may be approximated by a con- 
tact potential, there is a rapidly growing interest in a 
novel type of cold gases in which electric or magnetic 
dipole-dipole interactions (DDI) are crucial for the oc- 
curing phenomena. These so-called dipolar gases in- 
clude atoms with large magnetic moments [1-3], polar 
molecules [4-6] and Rydberg gases [7]. The distinct na- 
ture of the dipolar interactions leads to a wealth of novel 
physics, including a geometry-dependent stability and a 
roton-like minimum in the excitation spectrum (for re- 
views see e.g. [8, 9]). 

Interestingly, the long-range dipolar interactions result 
in an inherent non-local nature of dipolar gases, partic- 
ularly striking in deep optical lattices. For non-dipolar 
systems, gases trapped in different sites of a deep lattice 
do not interact with each other. Hence, for a vanish- 
ing inter-site hopping, different sites may be considered 
as independent, uncorrelated experiments. In contrast, 
inter-site DDI play a substantial role even in the absence 
of hopping. Recent lattice experiments have shown that 
the inter-site dipolar interactions are a key element in the 
dynamics [10], as well as in the stability and collapse of 
dipolar condensates [11, 12]. 

Faraday patterns constitute a paradigmatic example of 
pattern formation in periodically driven systems [13, 14] 
ranging from classical fluids [15], through multimode 
lasers [16] and superfluid Helium [17]. Interestingly, Fara- 
day patterns may be observed in BECs by modulating 
the nonlinearity resulting from the interatomic interac- 
tions [18-23], as shown in recent experiments [24]. Fara- 
day patterns in Bose-Einstein condensates (BECs) may 
be directly linked to the spectrum of elementary exci- 
tations, and in this sense provide an excellent insight 



into fundamental properties of the condensates. In non- 
dipolar gases the Faraday pattern selection is determined 
uniquely for each modulation frequency due to the mono- 
tonically growing character of the excitation energy [25] . 
Interestingly, this is not anymore the case for dipolar 
BECs with a roton-like minimum in the excitation spec- 
trum [26]. As a result, it has been shown that Fara- 
day patterns in two-dimensional (2D) dipolar conden- 
sates present remarkable qualitative novel features [27] 

In this paper, we analyze quasi-one-dimensional 
(quasi- ID) dipolar condensates with periodically driven 
dipolar interactions. We demonstrate that Faraday pat- 
terns provide a clear example of the non-trivial role 
of confinement dimensionality in dipolar gases, showing 
that in the presence of a roton-like minimum in the ex- 
citation spectrum, Faraday patterns in a quasi-lD trap 
differ significantly with respect to the 2D case [27]. More- 
over, Faraday patterns provide as well an excellent tool 
for a study of non-local effects in dipolar condensates, 
as we illustrate with two parallel quasi- ID BECs, in the 
absence of tunneling. The non-local dipolar interactions 
between both BECs lead to an unfolding of the excita- 
tion spectrum into symmetric and anti-symmetric modes 
with respect to the transposition of the two condensates. 
We show that, as a consequence, at a critical driving fre- 
quency a transition between correlated (symmetric) and 
anti-correlated (anti-symmetric) Faraday patterns in the 
two BECs occurs. For the critical driving the emergent 
Faraday pattern differs from one realization to another, 
resulting from a spontaneous symmetry breaking mech- 
anism. 

The paper is structured as follows. In Sec. II we intro- 
duce the model for periodically driven quasi- ID dipolar 
condensates. Section HI is devoted to Faraday patterns 
in a single quasi- ID BEC, with a focus on the differences 
as compared to 2D condensates. Section IV is dedicated 
to the effects of the inter-condensate dipolar interactions 
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on the Faraday pattern selection in two parallel disjoint 
quasi-lD dipolar condensates. We conclude in Sec. V. 



II. MODEL 

We consider in this paper quasi- ID dipolar BECs, 
either in a single trap (Sec. Ill) or in two parallel 
traps (Sec. IV). Since the former case may be consid- 
ered as a particular realization of the latter, we present 
in this Section the general formalism for parallel quasi- 
lD BECs aligned along the z axis, and separated along 
the y axis by a distance A. We assume the potential bar- 
rier separating both quasi- ID BECs sufficiently large to 
suppress any hopping between them. Each condensate 
experiences a strong harmonic confinement of frequency 
cj^ in the x-y plane and no confinement along the z di- 
rection. The atoms possess a magnetic dipole moment 
/i (the results are equally valid for electric dipoles) ori- 
ented by an external field along the y axis. We employ 
dimensionless expressions, using units of frequency uj^ 
and length = y/h/Muj^, with M the particle mass. 

Due to the strong x-y confinement we assume that the 
system remains in the ground state of the x-y harmonic 
oscillator (this condition is self-consistently verified), and 
employ the non-local non-linear Schrodinger formalism 
developed in Refs. [28, 29] for a stack of quasi-lD dipo- 
lar BECs, to obtain the coupled equations for the wave 
functions ^l^j{z) in traps j = 1, 2: 



e(q) 



idfi/jj (z) 



27r 



1 



dz^grijiz) 



^j{z). (1) 



Short-range interactions are characterized by the cou- 
pling constant g = g^^no/27rhuj^l^, where no is the linear 
density, and g^^ = Anasc^^ /M, with age the s-wave scat- 
tering length. The DDI are determined by the coupling 
constant gd = g^^nQ/27rhuj^l^^ where g^^ = /ioM^/47r, 
with /io the vacuum permeability. In Eq. (1), fimi^z) 
is the Fourier transform of the linear density n^(z) = 
\^l)m (^)l^, and 



dk 



ke i''^ 



[{k^ - 2kl) Jo {kAp) - 3fc2 Ja (kAp)] , 



(2) 



where Jn{x) are the Bessel functions of the ffist kind. 

In the following we consider a parametric modulation 
of the dipole-dipole interactions 



gd{t) =^^(l + 2acos(2a;t)). 



(3) 



where a characterizes the modulation strength. Such 
modulation may be implemented with intensity oscilla- 
tions of the polarizing electric field for the case of polar 
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FIG. 1. (Color online) Excitation spectrum e(^) of a single 
quasi-lD dipolar EEC, with g = -0.1007, gd = 0.0629. Note 
the roton minimum at cjr = e(^r) and the maxon maximum 
at (jJm- For a driving frequency uJr < oj < uJm there are 
three possible momentta ^1,2,3 obeying the resonance condi- 
tion e(^) = cj. 



molecules, or with additional transverse magnetic fields, 
which lead to a precession of the dipole moment orienta- 
tion, for the case of magnetic dipoles. 

The modulation of gd induces Faraday waves. With 
the aim of examining the growth of such patterns, we 
introduce the following ansatz for the wave functions: 



iljj {z,t) = iljj {l^Aj{t) cos{qz)) , 



(4) 



which describes well the physics of the pattern in the lin- 
ear regime, where the modulation is weak and we may 
consider each momentum component q of the pattern 
separately. In Eq. (4), we introduce the complex am- 
plitude Aj{t) = Uj{t) -\-ivj{t)^ which determines the per- 
turbation from the initial homogeneous solution tpj^ = 
ex.-p{—iJIj[t-\-{Qj/uj) sin(2a;t)]}, where = a{l—gd/'pj), 
and JIj = g-\- ^gd ^|m-j| (0) is the chemical potential. 
Inserting Eqs. (3) and (4), into Eq. (1), and linearizing 
in Aj we arrive at the system of equations describing the 
modulation dynamics 
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2g Uj 



47r >r^ 



= 0. 



(5) 



III. 



FARADAY PATTERNS IN A SINGLE 
QUASI-ID DIPOLAR EEC 



We consider in this Section the case of a single con- 
densate, being particularly interested in the differences 
between the emergent Faraday patterns in a quasi- ID 
trap and those predicted in Ref. [27] for a 2D condensate. 
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Employing a similar Bogoliubov analysis as the one pre- 
sented in Ref. [28] , we obtain the spectrum of elementary 
excitations in the considered case (see Fig. 1): 



e('/) = ^^(^ + 25+y5,Fo(g)) (6) 

where 

J^o(9) = l + ^9V^/'EiHV2), (7) 

with Ei{x) the exponential integral function. Using 
Eqs. (5) and (6), we arrive at the corresponding Mathieu 
equation [30]: 

+ [e^(^) + 2ujH{q, uj, a) cos(2cjt)] = (8) 

with 

h{q, cj, a) = ^g^aq'^Foiq). (9) 

Following Floquet theorem [31], the solutions of Eq. (8) 
are of the form u{t) = e^^f{t) where f{t) = f{t + tt/uj) 
and a{q^uj^a) is the Floquet characteristic exponent, 
which can be found numerically. If the real part a = 
Re(a) > 0, the homogeneous quasi-lD BEC becomes dy- 
namically unstable and Faraday patterns emerge. The 
typical wave length of the pattern will be determined by 
the most unstable mode, i.e., that with the largest a. In 
the limit of small driving amplitude, a ^ 0, the proper- 
ties of the pattern are governed by momenta q obeying 
parametric resonances en{q) = nuj. 

Contrary to non-dipolar BECs with a monotonic spec- 
trum e{q)^ dipolar gases may offer a more complex roton- 
maxon spectrum [26] (Fig. 1). As a consequence of this 
non- monotonic character, for a specific range of cj, be- 
tween the roton and maxon frequencies {uOr and uOrm re- 
spectively) there are three values qi < q2 < qs satisfying 
the resonance condition e{q) = uj. Figure 2 shows the 
stability diagram for a driving frequency in this partic- 
ular window. As expected, for small amplitudes a, the 
three instability tongues (white regions) correspond ex- 
actly to ^^1,2,3 (Fig. 1). This opens an interesting question 
about which of the three modes dominates the pattern 
formation. For a 2D geometry, Ref. [27] showed that 
when modulating dipole-dipole interactions, the most un- 
stable mode corresponds to the intermediate momentum 
Q2 < Qr^ with qr the roton momentum. Crucially, as 
we show below, this is not the case in a quasi-lD dipolar 
condensate. This striking contrast between quasi- ID and 
2D predictions illustrates once more the key role played 
by the trapping geometry in dipolar gases. 

The problem of the most unstable mode is best un- 
derstood employing a series expansion of the Floquet 
exponent with small parameter b{q^uj^a) [20, 32, 33], 
which, for the first parametric resonance e{q) = uj^ yields 
a :^ h{q^uj^a) /2 oc q^Fo^q). Remarkably, in contrast to 
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FIG. 2. (Color online) Stability diagram for the parameters 
of Fig. 1, as a function of the perturbation strength a and 
momentum q. The unstable region is depicted in white. The 
solid line indicates the most unstable mode, and the dashed 
line the roton momentum. 
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FIG. 3. (Color online) Most unstable momentum as a func- 
tion of the driving frequency uj for the parameters of Fig. 1, 
with a — 0.01. The horizontal dashed lines indicate the roton 
and maxon frequencies (cjr.m), and the vertical line the roton 
momentum. 



the 2D case, we find that the most unstable mode cor- 
responds to the largest momentum q^, > qr (solid line in 
Fig. 2). Fig. 3 depicts a momentum of the most unsta- 
ble mode as a function of the driving frequency uj. The 
plot confirms that for all frequencies within the window 
uJm < UJ < uJr the momentum characterizing the most 
unstable mode is larger than the roton momentum, con- 
tradicting the prediction for a 2D pancake geometry [27]. 
For UJ < ujr^ alike in the 2D case, the observed modula- 
tions are dominated by higher resonances with q in the 
vicinity of qr. However, unlike the 2D scenario, even in 
this regime the most unstable mode in a quasi- ID BEC is 
characterized hy q> qr. We emphasize that the different 



4 



nature of the Faraday pattern reported here stems solely 
from the quasi- ID character of the condensate, which 
leads to a specific momentum dependence of b{q^uj^a) 
that differs from that in 2D. 

We have simulated numerically the time evolution 
of the non-local non-linear Schrodinger equation (1) 
with the parametrically driven nonlinearity, according to 
Eq. (3). The emergent pattern has been examined by 
means of Fourier transform of the condensate density, 
which confirmed the results for the most unstable mode 
that we obtained within the Mathieu analysis. 



IV. FARADAY PATTERNS IN TWO ID 
DIPOLAR BECs 

We now turn to the study of Faraday patterns in two 
parallel quasi- ID dipolar BECs. For non-dipolar BECs, 
in the absence of hopping, each BEC behaves indepen- 
dently, and hence an experiment with two BECs reduces 
to two uncorrelated experiments with a single conden- 
sate. The situation is radically different in dipolar BECs, 
since, despite the absence of hopping, the non-local char- 
acter of the dipolar potential gives rise to a coupling 
between the two BECs, with the strength of the inter- 
condensate interactions governed by Fi{kz). These non- 
local interactions lead to a collective character of the 
elementary excitations that are shared among the two 
quasi-lD condensates [28, 34]. Consequently, the excita- 
tion spectrum unfolds into two branches 




gd{Fo{q)±F,{q))], (10) 



which correspond, respectively, to symmetric and anti- 
symmetric states with respect to the transposition of 
traps j = 1 ^ j = 2. 

Interestingly, this implies that a periodic modulation 
of the dipolar interactions yields two different parametric 
resonances for each driving frequency uj = €±{q±)^ even 
in the absence of the roton minimum (see Fig. 4). Note, 
that the patterns are characterized not only by their mo- 
mentum q± but also by their symmetric (+) or anti- 
symmetric (— ) character. In analogy to Sec. Ill, the dou- 
ble solution raises a fundamental question about which of 
these two modes is the most unstable, and hence provides 
the dominant Faraday pattern. We stress that this non- 
trivial physics stems directly from the inter-condensate 
interactions, which lead to the splitting between the two 
branches in the spectrum, being a qualitatively new fea- 
ture of dipolar condensates. 

Similarly to the previous section, we employ Eqs. (5) 
for j = 1,2, and the spectra (10). In turn, we obtain 
two decoupled Mathieu equations for the symmetric and 
anti-symmetric combinations u± = Ui ^ U2'' 



d?u± 



[e\{q) ^2uj%±{q,uj,a)co^{2ujt)\ u = {) (11) 




FIG. 4. (Color online) Elementary excitations of two paral- 
lel quasi-lD dipolar BECs, for g = -0.0629, = 0.1749, 
and A//^ = 6. Note the two branches of the elementary 
excitations e±(^), corresponding, respectively, to symmetric 
and anti-symmetric modes with respect to the transposition 
of traps j = 1 o j = 2. 



with 



h±{q,uj,a) 



27r 



g^aq^ {Fo{q) ± F,{q)) , 



(12) 



to which we apply the Floquet analysis employed in the 
study of Eq. (8). As in the case of a single BEC, the 
first parametric resonances e±{q±) = u are character- 
ized by the Floquet exponent a± o:^ b±{q,uj^ a) /2 oc 
q'^{Fo{q) ± Fi{q)), and the emerging Faraday pattern is 
determined, for each driving frequency separately, by the 
mode with the largest a. Remarkably, the involved mo- 
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FIG. 5. (Color online) Real part a± of the Floquet exponent, 
corresponding to the the first parametric resonance for the 
symmetric and anti-symmetric excitation branches e±{q) — 
cj, as a function of the driving frequency uj. Note that at 
a critical frequency ujc = 0.055, both exponents are equal, 
cr+ = cr_ , indicating a transition between the symmetric and 
the anti-symmetric Faraday pattern. In the figure we employ 
g = -0.0435, g^ = 0.0437, A = 6/±, and a = 0.02. 
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FIG. 6. (Color online) Analysis of the pattern selection as 
a function of the driving frequency cj, in the neighborhood 
of the critical frequency ujc (for the same parameters as in 
Fig. 5). The solid lines represent the excitation branches. For 
each uj we indicate with a circle a momentum value where 
the numerical Fourier transform nj{kz) of the Faraday pat- 
tern shows a clear maximum. For uj well below (above) uJc we 
observe a single peak at ^+ (q-) indicating that a symmet- 
ric (anti-symmetric) Faraday pattern emerges (see insets). In 
the vicinity of Uc (shaded region) , both modes are equally un- 
stable and we observe the two corresponding peaks occuring 
simultaneously in the Fourier transform (see text). 



mentum dependence of Fo(^) i^i {q) leads to an intricate 
relation between the Floquet exponents and the driving 
frequency cj, as presented in Fig. 5. 

Crucially, the curves (J±{oj) cross at a critical fre- 
quency UJc- In consequence, we expect a distinct tran- 
sition, as a function of the driving frequency cj, be- 
tween the symmetric Faraday pattern for u < uJc 
and the anti-symmetric pattern for uj > uj^ Such 
transition is marked by an abrupt change of the pat- 
terns from a maximum- maximum alignement (correlated 
patterns) to a maximum- minimum alignement (anti- 
correlated patterns), as depicted in the corresponding 
insets of Figs. 6 and 7. Moreover, for uj = uJc the pat- 
terns in both condensates exhibit a pronounced change 
of the wavelength of the modulation, from l/qj^{ujc) to 

{uJc). 

This transition has been confirmed by means of direct 
numerical simulations of Eqs. (1), with the parametric 
driving governed by Eq. (3). As for a single condensate, 
we Fourier transform the density of each condensate to 
obtain the dominant momenta of the emergent Faraday 
patterns. The results, in the vicinity of the critical fre- 
quency cjc, are depicted in Fig. 6, where, on top of the 
spectra e± , for each driving frequency uj we indicate with 
a circle the momentum value where the numerically eval- 
uated fijikz) shows a marked maximum. We find that, in 
agreement with the results for (t±{uj) presented in Fig. 5, 
for UJ well below uJc the pattern presents a single momen- 
tum component at qj^^ being characterized by a corre- 



lation between the patterns in both quasi- ID BECs. In 
contrast, for uj well above uJc a single momentum compo- 
nent q- is observed, and the patterns in the two quasi- ID 
BECs are anti-correlated. 

In order to quantify the transition between correlated 
and anti-correlated patterns we introduce the correlation 
coefficient 



^dzSnAz)-Sn,{z) 



JdzSl{z)-JjdzSl{z) 



(13) 



where Snj {z) = rij {z) — nj, with nj the average density 
in a trap j. Pattern correlation is then characterized by 
r > 0, whereas anti-correlation leads to r < 0. Fig. 7 il- 
lustrates the radically different time evolution of the cor- 
relation coefficient below and above the critical driving 
UJc. Clearly, for frequencies sufficiently smaller (larger) 
than UJc the system arrives at perfectly correlated (anti- 
correlated) pattern with |r| = 1. 

An interesting scenario occurs for driving frequencies 
in the vicinity of the critical uJc (shaded region in Figs. 5 
and 6), where both the symmetric pattern with wave- 
length l/q_^{uJc) and the anti- symmetric pattern with 
wavelength l/q-{uJc) are equally unstable. As a result, 
the Fourier transform of the density in each quasi- ID 
BEC shows a simultaneous appearance of both momen- 
tum peaks, g+ and q- (see Fig. 6). 

Note that at = cjc, not only e+(g'+) = e_(g'_) but 
also 6+ (g+ , cj, a) = b- {q- , cj, a) and hence the two Math- 
ieu equations (11) for U-^ and U- become identical. This 
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FIG. 7. (Color online) Correlation function r(t) for the same 
parameters as in Fig. 6 {uJc = 0.055). The upper curve, which 
corresponds to cj = 0.025 < cjc, approaches r = 1 indicating 
a perfectly correlated pattern in both quasi- ID traps. The 
lower curve, which corresponds to cj = 0.08 > cjc, reaches 
r — —1 proving a perfect anti-correlation between the Fara- 
day patterns in the two traps. The insets show the corre- 
sponding numerical results for the density distribution nj{z) 
for t = 11000, with the bright (dark) colors indicating den- 
sity maxima (minima). Naturally, for sufficiently long times, 
well beyond the linear regime, the Faraday patterns and their 
correlations are eventually destroyed. 
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FIG. 8. (Color online) Population imbalance x(t) between 
the two peaks at q±{(jOc) for the critical driving a; = cjc (for 
the same parameters as Fig. 6). Note the 2uj periodicity 
(T = tt/cj = 57.1) that stems from a spontaneous symme- 
try breaking mechanism (see text). 



symmetry is however spontaneously broken in experi- 
ments due to quantum and thermal fluctuations, which 
lead to different initial conditions (populations) for both 
modes, that change randomly from one realization to 
another. This spontaneous symmetry breaking mecha- 
nism is best studied quantitatively by considering the 
relative weight of the momentum peaks at (7+ and q- in 
the Fourier-transform of the density h(kz). To this end, 
we define the imbalance parameter 



x{t) 



h{q+,t) - n{q-,t) 
h{q+,t) +h{q-,t)' 



(14) 



For CO well below or above cjc, once the pattern emerges, 
x{t) = ±1. In the vicinity of cJc, however, the imbal- 
ance parameter x(t) shows a clear periodicity with fre- 
quency 2uj (see Fig. 8). Note that these oscillations do 
not result from non-linear competition, as they occur 
well within the linear regime. In fact, the 2uj oscilla- 



tions of x{t) originate in different, spontanously chosen, 
initial conditions for U-^ and which lead to their dif- 
ferent time evolution that can be well approximated by 
u±{t) = {u^cos{ujt) + sin(co't)) exp((jCo't), where 
are the constants determined by the initial conditions. 
Furthermore, spontaneous symmetry-breaking leads to a 
different result for the imbalance x{t) from one realiza- 
tion to another, what we have confirmed by considering 
small random differences in the initial conditions for our 
numerical simulations of Eqs. (1). 



V. CONCLUSIONS 

Faraday patterns in dipolar BECs are crucially depen- 
dent on the unique properties of the dipole-dipole inter- 
actions. In particular, due to the long-range anisotropic 
nature of the dipolar interactions, the character of the 
Faraday patterns depends strongly on the dimensional- 
ity of the condensates. We have shown that for period- 
ically modulated dipolar interactions, Faraday patterns 
in 2D and ID geometries differ substantially in the pres- 
ence of a roton minimum in the excitation spectrum. 
Moreover, for parallel quasi- ID dipolar BECs, the inter- 
condensate interactions lead, even in the absence of hop- 
ping, to an excitation spectrum characterized by sym- 
metric and anti-symmetric modes. This, in turn, gives 
rise at a critical driving frequency to a marked transition 
between correlated and anti-correlated Faraday patterns 
in the two condensates. Interestingly, at this transition 
point the Faraday pattern selection stems from a spon- 
taneous symmetry breaking mechanism. 
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